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ABSTRACT 

In this paper, we examine in detail the key structural properties of high redshift dark 
matter haloes as a function of their spin parameter. We perform and analyze high res- 
olution cosmological simulations of the formation of structure in a LCDM Universe. 
We study the mass function, shapes, density profiles, and rotation curves for a large 
sample of dark matter haloes from z = 15 — 6. We also present detailed convergence 
tests for individual haloes. We find that high spin haloes have stronger clustering 
strengths (up to 25%) at all mass and redshift ranges at these early epochs. High 
redshift spherical haloes are also up to 50% more clustered than extremely aspherical 
haloes. High spin haloes at these redshifts are also preferentially found in high density 
environments, and have more neighbors than their low spin counterparts. We report 
a systematic offset in the peak of the circular velocity curves for high and low spin 
haloes of the same mass. Therefore, estimating halo masses without knowledge of the 
spin, using only the circular velocity can yield errors of up to 40%. The significant de- 
pendence of key structural properties on spin that we report here likely has important 
implications for studies of star formation and feedback from these galaxies. 

Key words: cosmology: dark matter - cosmology: early Universe ~ galaxies: high- 
redshift - galaxies: formation 



1 INTRODUCTION 

The currently favored model that describes the formation of 
structure in the Universe is the A cold dark matter (LCDM) 
paradigm. In this model, the initial density distribution 
of the Universe was nearly homogenous, with small Gaus- 
sian density perturbations imprinted during an inflationary 
epoch. These fluctuations expand linearly, until the over- 
dense regions undergo non-linear gravitational collapse to 
form bound dark matter haloes. These haloes form in a hi- 
erarchical fashion: small haloes form first, and then larger 
ones assemble later via merging. In the LCDM paradigm, 
baryons follow the dark matter. Since they can dissipate 
and cool, baryons condense, and eventually form observable 
galaxies in the centres of dark matter haloes. 

The properties of dark matter haloes in the context of 
the LCDM paradigm have been studied in detail using nu- 



growth of structure in the Universe. These theoretical stud- 
ies provide the framework within which the role of baryons 
and details of galaxy formation can be probed. While coUi- 
sionless dark matter in the LCDM paradigm interacts only 
gravitationally, baryons dissipate, have pressure, cool, form 
stars, and interact with radiation. These, and other ef- 
fects, introduce complications when trying to understand 
the properties of dark matter haloes such as their mass, 
angular momentum, shape, and density profiles from ob- 
servations of the baryonic component. There are, however 
two techniques that have allowed a more direct probe of 
the dark matter: gravitational lensing obser vations (e.g . 
Fischer et allbOOOl: [McKay et al...2002.:.Hoekstra et al.ll2003 



Mandelbaum et al.ll2006yLimousin et al.ll2007l: IParker et all 



20071 : I Evans fc Bridle! l2009ll. and measurements of galax; 



rotation curves 



uuai. and rneasure r nents ot galaxy 
Rubin etal] Il985l: iTrimbld IiQStI : 



creasing resolution (e.g. 


Davis et al.lll985; Frenk et al. 1988 


Efstathiou et al. 


1988; 


Katz, Hernauist & Weinberg 


199C 


KaufFmann et al 


.. .19991: 


iBuUock et al.l |2001B|: iFrenk 


2002 



Springcl 200i). This approach has been very fruitful in 
providing us with a detailed picture of the assembly and 



IPersic. Salucci. fc St el"l996':' Salucci et al.H 2007'). Due to the 
difficulties and assumptions required to translate the ob- 
served baryonic properties to dark matter halo properties, 
cosmological N-body simulations offer a powerful tool to un- 
derstand the properties and statistics of the dark matter 
haloes. 

Even with dark matter only numerical simulations, 



2 Davis & Natarajan 



much has been learned about the assembly of dark mat- 
ter haloes, including the halo mass function, halo clustering, 



^^^^^ ^^^^1 ^^^^ ^ ^^^^ ^^^.^^^^^^ -"J-l-l • . ^^^^^^ 

2002; Reed et al.[2003. 2007; Bctt et al. 2007; Gao & Whita 


2007 




Reed et alj 2009; 


Maccio, Button. & van den BoschI 


200S 


Faltenbacher & Whites i2010l). However, there have 



been few detailed studies of dark matter halo properties at 
high redshifts. This is partly due to the number of particles 
required to resolve high redshift, low mass haloes, and still 
match observations of larger haloes at lower redshifts. These 
restrictions until recently prevented the detailed study of a 
statistically significant sample of collapsed haloes at high 
redshifts. As the observational frontier is pushed to higher 
and higher redsh ifts with reports o f the detection of galaxies 
out to z ~ 7—8 l|Oesch et al.ll2O10l ). a deeper understanding 
of the properties of the dark matter haloes that host these 
most distant galaxies is critical as well as extremely timely. 

A few recent studies have exam ined specific dark matter 
halo propertie s at h igher r edshifts. iHeitmann et al] (|2006l ). 
I Warren et all l|2006l ). and iReed et al.l (|2007l) focus on the 



mass function of high redshift haloes. iMoore et al.l (|2006| ) 
trace the spatial distribution of dark matter halos from 
z = 12 to the present d ay to understand their effect on 
galaxy mass haloes today. IJang-Condell fc Hernguistl (|200ll ) 
use low resolution simulations to determine the spin and 
shape para meters of dark m atter haloes at 2; = 10. In a re- 
cent study (jPavis fc Natarajan .20091) we reported the re- 
sults of the first high redshift and high resolution study 
to follow the growth of angular momentum in dark mat- 
ter haloes in the mass range 10^ M© to 10* Mq from z — 15 
to z = 6, a period spanning 700 Myrs of cosmic time. We 
found that the spin distribution at these early epochs can 
be fit by a log-normal distribution as at lower redshifts. In 
addition, we examined the two-point correlation function of 
haloes and found that higher spin haloes are more clustered 
by factors up to 25% compared to their low spin counterparts 
at a given mass. This finding extended across all mass bins 
and redshifts in our previous study, i.e. from lO*" — 10* Mq 
and from z — 15 — 6. 

This paper builds on our earlier work by investi- 
gating the role angular momentum and the environment 
play in the determination of structural properties of 
dark matter haloes at these epochs. In the LCDM 
paradigm, haloes acquire an gular mome n tum by tidal 
torques from their neighbo rs iHovld I1949I ; iPeebled 1 19691 ; 
iDoroshkeviciil I197G| ; IWhitd 1 19841 ) This picture for the 
acquisition and growth of angular momentum has been 
shown to be generally accurate in N-body simulations 
wherei n angular momentum init ially grows linearly with 
time (|Barnes fc Efstathiou 1987 ) and then slows down 



at l ater times ( Sugerman. Summers, fc Kamionkowskil 
l2000h . Linear theory, however, overpredicts the an- 
gular momentum when compared to the fully non- 
linear N-body simulations ([Barnes fc Efstathioul 



19871; ISugerman. Summers, fc Kamionkowskil 2000l ; 



Porciani. Dekel fc Hoffman 
Vitvitska et all (|2002D 



2OO2I ). 



point out. 



In addition, 
linear theory 



as 
pre- 



dicts the angular momentum of a halo at a given redshift, 
but not the angular momentum of any particular progenitor 
at an earlier redshift. Thus, it becomes impossible with 
linear theory to trace the evolution of a halo's angular 
momentum in a hierarchical Universe evolving via mergers. 



Vitvitska et all 1I2OO2I) ; i Maller. Dekel. fc Somer^dii^ l|2002h ; 
Hetznecker fc BurkertI ( 20061 ) all note that mergers do 



affect the sp in of the halo in add ition to the tidal torque 
model. iD'Qnghia fc Navarrol \200i ) study mergers and spin 
evolution explicitly and argue that mergers only affect the 
spin of unrelaxed haloes, and find that relaxed, isolated 
haloes show no correlation between spin and merger history. 

One way to study the acquisition of angular momen- 
tum is to correlate information about the environment with 
halo properties. Previous studies have shown that halo clus- 
tering strength depen ds on the a ngular momentum of the 
halo at low redshifts (iBett et al.ir2 007; Gao & Whit3 j2007l ; 
iFaltenbacher fc Whitell2010l ). lAvila- Reese et al.i (2005) find 
that galaxy mass haloes inside clus ters have smaller spins 
than haloes in the field or voids, and lReed et al.l l|2005h find 
low specific angular momentum in subhaloes near the cen- 
tral host halo. 

Observations using the Sloan Digital Sky Survey 
(SDSS) show that the spin parameter, 



|l/2 



(1) 



where J is the total angular momentum, Etot is the to- 
tal energy, and M the halo mass, has little dependence 
on th e local environmental density (|Cervantes-Sodi et al.l 
I2OO8I). It remains unclear, however, whether the spin pa- 
rameter derived from the baryonic disk model used in in- 
terpreting SDSS data correlates well with the host dark 
matter halo's spin which is what is assumed. These results 
are found using galaxy neighbors t o trace the lar ge scale 
tidal field ( Cerv antes-Sodi. Hernandez fc Parkll2O10l ). As an 
example of the difficulties in relating baryonic properties 
to the host dark matter halo, iQuadri. Moller. fc NataraianI 
(j2003l ) report how slight spatial offsets between the dark 
matter and the baryonic disk create disturbed lensing con- 
figurations which can be easily misinterpreted if the mis- 
alignment is not in c luded in the mass reconstruction. Also, 
iMandelbaum et al.l (|2006l ) report the misalignment of light 
ellipticity with halo ellipticity in the SDSS catalog. These 
findings illustrate the complexities when inferring dark mat- 
ter halo properties from baryonic observations. 

Assembly bias refers to the observation that the cluster- 
ing strength of dark matter haloes depends on an additional 
parameter beyond just halo mass. Assembly bias has been 
studied in simulations by examining halo formation time, 
concentration, sha pe, triaxiality, veloc i ty structure, and sub- 
structure content jHarke r etal? ' 20(361; IWechsler e t aP '200d; 
Bett et al.ll2007l;lGao fc W hite 20 071; I Jing, Suto fc Mo 200j; 



Wetzel et all |2007| ; I An'gu lo. Baug h. fc Lacevll2008i ). These 



previous works show convincingly that halo clustering de- 
pends on more than just halo mass. However, all of these 
studies have been at low redshifts {z < 5) and for massive 
haloes (M > 1 x 10^° Mq). 

In this paper, we extend previous work by examining 
the formation and growth of dark matter haloes at high red- 
shift, with an emphasis on studying the role angular momen- 
tum and environment play in regulating the structural prop- 
erties of dark matter haloes. We limit ourselves to haloes in 
the mass range of lO'' M© to a few times 10^ M© , and the 
redshift range z = 15 to z — 6. This allows us to focus on the 
dark matter haloes that will likely host the first generation 
of stars and galaxies. 
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We outline our paper as follows. In Section [2] we sum- 
marize studies of angular momentum and assembly bias at 
low redshift to provide the framework for our findings at 
high redshift. We describe our simulations in Section |3] and 
present the results of convergence tests in Section |4l Our re- 
sults from the correlation of the spin parameter to the halo 
environment are presented in Section O and that of the ef- 
fect of angular momentum on halo structure in Section [G] 
We conclude with a discussion of the implications of our 
results for high redshift galaxy formation. 



2 PREVIOUS STUDIES AT LOW REDSHIFT 

In this section, we summarize earlier findings pertaining to 
measurements of halo spin and clustering to provide the 
context for our findings. Many numerical simulations have 
shown that for massive haloes at low redshift, the distri- 
bution of the dimensionless spin parameter follows a log- 
normal distribution, 



P(A) = 



Av/27ro- 



exp 



-(ln(A/Ao))^ 
2a2 



with tjrpical val ues of Ao 
(e.g., rP'Onghi a fc Navarro 



(2) 



0.5 



Bailin fc Steinmetzl 



2007 



199e 



^ 03 5 and a 

I2OO7I: iBett et all 
■20051: ICole fc LacevI 
Steinmetz fc Bartelmannlll995l : [Warren et al.lll992l '). 

Previous studies have also shown that dark matter 
haloes are ge nerally t riaxial w ith a preference fo r pro- 



lateness (e.g. iBett et al . 2007; iBailiii fc Steinmetzl l2005l: 



Faltenbacher et al.ll20oi ICole fc LacevI 1 19961 : 1 Warren et all 



1992l : lFrenk et al.Hl988f ). iBett et all (|2007l ') find that nearly 



spherical haloes have smaller spins, while there is only 
a weak trend of halo triaxiality with spin. In studying 
halo concentration versus spin, it is seen that when un- 
relaxed haloes are eliminated from the sample there is 
only a very weak sur viving (if any) co r relation between 
these two parameter s (|Macci6 et al.ll2007l : iNeto et al.ll2007l : 
iBuUock et JlbOOlBl) . 

iBett et all ||2007| ) : iGao fc Whitd (|2007l ) and 
iFaltenbacher fc Whitd l|201oi r all find that haloes with 
larger spins are more clus tered than low sp i n halo es at 
a given mass. H owever, lAvila-Reese et all (|2005l ) and 
iReed et al.l (|2005') find in their simulations that haloes in 
cluster environments have smaller spins and are more spher- 
ical th an their counterparts in the field, and iHahn et al.l 
i2003) find that haloes in filaments have larg er spins than 
haloes of the same mass in clusters or in voids. iMaccio et al.l 
l|2007h find that there is no environmental dependence on 
the spin parameter for haloes at a given mass. Thus, there 
appears to be some question as to the extent to which the 
environment affects the angular momentum properties of 
dark matter haloes. In this work, we explore this relation- 
ship at high redshift using both the clustering strength and 
the local density to characterize the environment. 



3 DESCRIPTION OF SIMULATIONS 

We run a series of N-body simulations to follow the growth of 
dark matter haloes from z ~ 100 down to z = 6. We choose 
the particle mass such that a 10^ Mq /h dark matter halo 



has 100 particles. For 512 particles, this requirement sets 
the comoving box size at 2.46Mpc/h and the particle mass 
at Mum ~ 1.0 x lO*M0/h. The initial condi tions are gen- 
erated using a parallelized version of Grafic (jPrunet et al.l 
[2008), which calculates the Gaussian ra ndom field for the 
dark matter particles. We use Gadget-2 (|Springelll2005l ') to 
follow dark matter particles down to a redshift of 2 = 6, 
with output snapshots at 2 = 15, 12, 11, 10, 9, 8, 7, 6. We use 
the WMAP3 (|»m, »a, Ob, ft, n, as j = 10 .238, 0.762, 0.0416, 
0.732, 0.958, 0.7611. ISpergel et al.ll2007l ) and the WMAP5 
({nM, nA,nb,/t,n,cr8} = 1 0.258, 0.742, 0.044, 0.719, 0.963, 
0.796]-. [Dunklev et al.ll2009l ) cosmological parameters for our 
simulations. The WMAP3 cosmology was used for our first 
runs studying numerical convergence of our measurements 
of the angular momentum, and the WMAP5 cosmology was 
used for the results presented in Sections [5] and [B] Table 
1 shows the runs and the relevant parameters used in this 
paper. Gadget-2 uses a softening length, e, to soften the 
gravitational force to prevent spurious 2-body interactions. 

To identify collapsed dark matter h aloes, we use the 
publicly available HOP code provided bv lEisenstein fc Hut! 
(1998). This method groups particles with their densest 
neighbor. After grouping, density thresholds are used to en- 
sure that haloes are not being over counted due to a halo 
being a subhalo within a larger overdensity. We choose the 
density thresholds in order to match th e high redshift mass 
function described in iReed et all (|2007| ). 

Once the haloes are identified, each particle in the halo 
is tested to s ee if i t is actually bound to the halo. We use 
SKID (Stadel 200lh to do the unbinding. SKID finds the po- 
tential and kinetic energies for all particles in a given halo 
and removes the most unbound particle from the halo. Suc- 
cessive iterations are performed until all particles are either 
bound or there are no more particles in the halo. Without 
unbinding, angular momentum properties could easily be 
dominated by a few transient particles not representative 
of the collapsed halo. We then define the halo mass as the 
total mass of all the particles assigned to the halo. Having 
calculated masses for our haloes, we can measure the mass 
function of our sample, and fi nd that the mass function of 
haloes does not match that of iReed et al.l (|2007l ) at the low 
mass end unless particles are unbound. We have in our halo 
catalogue ~ 24, 200 haloes in the mass range 10^^°' ^ Mq at 
z = 6, and « 16, 100 at z = 10. In the mass range 10''='=°'^ Mq 
we have ~ 2, 600 haloes at 2 = 6 and ~ 1250 at z = 10. 

To ensure that our halo sample is representative, we 
verify the halo mass function of our runs against theoreti- 
cal predictions, as we ll as against other simula tions. We use 
the Pres s-Schechter JPress fc Schechteilll974h and Sheth- 
Tormen (jSheth fc TormenI ll999F mass functions, and the 
fitting function provided bv lReed et ah! (|2007l ) derived from 
their simulations. We show in Figure [l]the results from the 
MedRes (which uses the WMAP 3 cosmology) run and the 
WM5 run. In both cosmologies, we find the mass function is 
poorly fit by the Press-Schechter mass fu nction. As the den - 
sity thresholds were chosen to match the lReed et al.l (|2007l ) 
mass function, it is unsurprising that it provides a better fit 
than the Press-Schechter function. However, we note that 
the Sheth-Tormen func tion provide s a go od fit: one that is 
slightly better than the lReed et al.l ([2003) function. 

In order to calculate A for our halos, we must first calcu- 
late the kinetic and potential energies as well as the angular 
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Name 




L 

[Mpc/h] 


e 

[kpc/h] 


A/dm 
[lO^Me/h] 


Cosmology 


LoRes 


256 


2.46 


0.3844 


8.0 


WMAP3 


MedRes 


512 


2.46 


0.18 


1.0 


WMAP3 


HiRes 


1024 


2.46 


0.095 


0.125 


WMAP3 


WM5 


512 


2.46 


0.09 


1.0 


WMAP5 



Table 1. Catalog of runs used in this paper. Np refers to the total number of particles, L to the comoving length of the box, e to 
the softening length, and Mdm to the dark matter particle mass. The assumed cosmological parameters are (Cm i f^Ai ^bi 'ii '''s) are: 
WMAP3 = (0.238, 0.762, 0.0416, 0.732, 0.958, 0.761); WMAP5 = (0.258, 0.742, 0.044, 0.719, 0.963, 0.796). We note that the HiRes run 
achieves the 1024'^ resolution only over one-eighth of the box, due to computational limitations. 
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56789 10 6789 10 
[!h] log M [Msun/h] log M [Msun/h] 

Figure 1. The hal o mass func t ion fo r the MedRes (right) and WM5 (left) runs plotted as diamonds. The solid (black) line is the 
mass function from iReed et al ] l|2007l 'l. the dotted (red) line is the Sheth-Tormen mass function, and the dashed (blue) line is the 
Press- Schechter mass function. 



momentum from the particles. The kinetic energy, Ek for a 
halo is the sum of the kinetic energy {^mv-v) of each par- 
ticle assigned to the halo. A halo's angular momentum, J, 
is calculated similarly as the sum of each particle's angular 
momentum (m fxv). The potential energy, Eg, of a halo 
is calculated using a direct summation: 

where G is Newton's gravitational constant, A'^^ the number 
of particles in the halo, and fij is the distance between par- 
ticles i and j. In all of these definitions, the velocities are 
with respect to the halo's mean velocity, and the positions 
are with respect to the centre of the halo. Throughout this 
paper, the centre of the halo refers to the location of the 
densest particle, which we use as a proxy for the location of 



the minimum of the halo's potential well. The only excep- 
tion will be in section [4l in which we use the centre of mass 
to cross reference halos in simulations of different resolution. 

Having now calculated Ek and Eq for our haloes, we 
can measure how virialized the haloes are at these early 
epochs. The scalar virial theorem states that for an iso- 
lated, coUisionless system in a steady state, the total kinetic 
energy should be equal to half the total potential energy: 
2Ek + Eg = 0. Thus, measuring the total kinetic and po- 
tential energies gives us insight into the dynamical state of 
these dark matter haloes. However, we note that the two 
key assumptions of the virial theorem (an isolated halo and 
steady state) are not strictly valid for these halos at these 
epochs. First, these haloes are still actively merging and ac- 
creting matter, and so have an effective surface pressure that 
adds to the total energy of the system. In addition to the 
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10' 
Halo Mass 



Figure 2. Ratio of 2E}^/\Eq\ for all haloes in the MedRes run at 
2 = 6 as a function of mass. For fully virialized haloes, this ratio 
should be 1. However, we find a mean value of 1.3, with only 0.2% 
of haloes having a value less than unity. The diamonds depict the 
mean ratio as a function of halo mass. The contours enclose ten, 
thirty, fifty, seventy, and ninety percent of the haloes (the same 
percentiles hold in all further figures with contours) 



surface pressure, we also note that the haloes are not nec- 
essarily in a steady state. The density of matter increases 
with redshift proportional to (l + z)^, and so at a redshift of 
z = 6, the mean matter density of the Universe is higher by 
a factor of almost 350. The higher density implies a larger 
merger rate for these dark matter haloes and hence a shorter 
time between mergers. Thus haloes are less likely to be in 
a steady state at high redshift. Finally, the virial theorem 
relates time-averaged values for the kinetic and potential en- 
ergies. However, we do not have the time resolution required 
for such a calculation. We therefore use instantaneous mea- 
surements of the kinetic and potential energies, with the 
understanding that this will induce scatter in the energies 
of our haloes. 

In Figure[2] we show the ratio of 2Ek/\Eg\ for all haloes 
in the MedRes run at z — 6. For a virialized dark matter 
halo, this should be approximately equal to 1. However, we 
find that the median ratio is 1.3, with fewer than 0.2% hav- 
ing a virial ratio less than unity at 2: = 6. At higher redshifts, 
the median ratio increases slightly to 1.5 a.t z — 10. We note 
that Hctzne ckcr & Burkert (200fS) found similar results, al- 
beit at lower redshifts. These results imply that very few of 
our haloes are actu ally virialized eval uated using the above 
definition. However, iBett et akl |20o3) la beled haloes with a 
ratio between 0.5 and 1.5 as relaxed, and iNeto et al.l (|2007l ') 
set the upper limit as 1.35. We find that 95% of o ur haloes fit 
thelBett et al.l (|2007l 'l criterion, while 73% fit the lNeto et~al] 
l|2007h criterion. 



4 CONVERGENCE AND NUMERICAL 
RESOLUTION EFFECTS 

In order to verify that we have the necessary particle resolu- 
tion to capture the angular momentum properties of haloes 



we ran three runs that had identical initial conditions, but 
different resolutions. This was achieved by creating the ini- 
tial conditions for a 1024^ sized run (HiRes) and then lump- 
ing particles to create a 512"* (MedRes) and 256'^ (LoRes) 
set of initial conditions. Thus, each particle in the MedRes 
run has the total mass of 8 particles in the HiRes run, and is 
assigned the average position and velocity of those 8 parti- 
cles. The LoRes run is similarly averaged out of the MedRes 
run. We note, however, that due to computational limita- 
tions, we were unable to run the entire box at the highest 
resolution. Therefore, we kept only one-eighth of the box at 
the highest resolution, and lumped the rest of the box to 
the medium resolution. In order to avoid biases due to the 
interaction of differing mass particles in the HiRes run, we 
only kept haloes which were entirely composed of the high 
resolution particles. This ensures that the haloes used for 
comparison are representative of the highest resolution. We 
used the centre of mass to cross-match haloes across the 
three runs. We match 85% of the haloes between runs. The 
unmatched haloes are likely cases where the group finder has 
joined neighboring groups with a tenuous bridge at a higher 
resolution that does not exist at lower resolution. In this sit- 
uation, the centre of mass would be very different between 
resolutions. 

In Figure [3] we compare the halo mass, spin, kinetic 
energy, potential energy, and total angular momentum be- 
tween the LoRes and MedRes simulations. The bottom right 
panel shows the cosine of the angle between the angular 
momentum vectors in the two simulations. Table[5]hsts the 
mean and standard deviation of the fractional difference of 
those same quantities as a function of number of particles 
in the LoRes run. We also include the offset in the centre of 
mass position, which was used to cross-match the haloes. We 
find good agreement in the masses of individual haloes, with 
a mean fractional difference of 14% at the lowest particle 
resolution (haloes with less than 1,000 particles). However, 
there is some spread in calculating the kinetic and potential 
energies, with differences up to 30% in the smallest bin. The 
angular momentum has the largest variation of the quan- 
tities used to calculate A. For the smallest haloes, we find 
a mean fractional difference of 73%, which only decreases 
to 25% in the largest haloes (haloes with more than 10, 000 
particles. We also find that the direction of the angular mo- 
mentum vector is biased at low particle resolution. 

These findings lead to a large spread in the spin param- 
eter for the same haloes at differing resolution. There does 
not appear to be any systematic offset in A between resolu- 
tions, as the top rig ht panel of Figure [3] shows. This is unlike 
I Trent! et al.l (|20ld ) who found that for haloes with less than 
100 particles, the spin parameter measurement is biased high 
at lower resolution. While we do not show the comparison 
between the MedRes and the HiRes subregion, we find the 
same trends with increasing resolution. One possible expla- 
nation for the larger dispersion in A may be the difficulty in 
defining an outer edge or boundary for dark matter haloes. 
A particle at the edge adds little to the halo's mass, but 
it will add a considerable amount of angular momentum. 
Thus, a measurement of any particular halo's angular mo- 
mentum or spin parameter will be dependent on the halo 
finding algorithm used in the study (jTrenti et al.ll2O10l ). 

In addition to the boundary issue, particle discrete- 
ness and resolution is also problematic when measuring halo 
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10^ 10^ 10^ 10^° 10" 10^^ 10^ 10^ 10^ 10^° 10" 10^^ 

E^-LoRes [Msun (km/s)^] Eg-LoRes [Msun (km/s)^] 




Jtot-LoRes [Msun kpc km/s] Mass-LoRes [Msun] 



Figure 3. Comparison of lialo properties in the LoRes and MedRes runs. The upper row plots haloes mass and spin, the middle row 
shows the kinetic and potential energies, and the bottom row shows the total angular momentum and the cosine of the angle between 
the angular momentum vectors in the two resolutions as a function of halo mass. Crosses (black) refer to haloes with fewer than 1000 
particles in the LoRes run, stars (blue) to haloes with fewer than 3160 particles, diamonds (red) to haloes with fewer than 10000 particles, 
and triangles (green) to haloes with more than 10000 particles. 
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Figure 4. Comparison of the spin distribution of all three reso- 
lution runs. The solid (black) curve shows the log-normal fit for 
the MedRes run, the dashed (blue) curve is the fit for the LoRes 
run, and the dot-dashed (red) curve is the fit for the HiRes run. 
Note that while the spin values for individual haloes have not 
converged, the statistical ensemble shows little variance between 
the three different resolutions. 



properties such as spin and angular momentum. As seen in 
Table [21 the dispersion decreases with increasing number 
of particles across all quantities. This implies that a large 
source of the scatter is due to particle resolution. 

The ensemble properties of the haloes' spin parame- 
ter reaches convergence at the medium resolution. Figure 
[4] shows the log-normal fit to the spin distribution for the 
three runs. The three runs have a mean A of 0.039 (LoRes), 
0.040 (MedRes), and 0.038 (HiRes). In the further analysis 
reported below, we use the results from the WM5 run, which 
has the same mass resolution as the MedRes run. 



5 SPIN AND ENVIRONMENT 



In our previous work l|Davis &: Nataraianl [20091 ). we found 
that at high redshift, higher spin haloes are more clustered 
than lower spin haloes at a given mass, using the halo-halo 
correlation function as a proxy for the environment. As in 
our previous work, we again define distinct spin bins from 
the lognormal distribution, each with roughly one third of 
the haloes, and denote them as high, median, and low spin 
bins. Using several different ways to characterize the en- 
vironment, we report our findings of correlations between 
measures of the environment with the spin parameter for 
haloes in the mass range 10'' M© to 10* Mq. 

We first use the distance to close neighbors as a measure 
of environment. This distance is calculated as the separation 
between each halo's centre (i.e. the densest particle). We 
measure both the distance to the 3'^'' (-D3) and T*** {D7) near- 
est neighbors. Both distances show the same result: haloes 
with closer neighbors have a small tendency to have larger 
spins. The second method that we use to measure environ- 
ment is the over-density within ten times the virial radius. 
We use as a proxy for the virial radius i?i78: the radius of 
a spherical region which has an average density of 178pcrit- 
The factor 178 is the overdensity criterion A{z) evaluated 



at z = 6. iBrvan fc NormanI (| 19981 ) give a fitting formula in 
a LCDM universe for A(2), but at high redshifts, A{z) is 
approximately the same as it would be in an Einstein-de 
Sitter universe, such that A{z) = IStt^ ~ 178. To calculate 
this radius, we use all particles surrounding the centre of the 
potential for the given halo, regardless of whether or not the 
particles are part of the halo. This is to account for the fact 
that not all haloes extended out far enough to reach the re- 
quired density threshold. In other words, some haloes have 
Ri7s outside the most distant particle assigned to the group. 
We then calculate the over-density within 107?i78, and find 
that haloes in over-dense regions have a slight tendency to 
have higher spins. Figure [5] shows the trends with spin for 
both measures of environment in two mass bins: io^='='''^ Mq 
and 10^^'''^ M0, and histograms for the lowest mass bin af- 
ter separating our sample by spin. The histograms for the 
10^ M0 and 10* M© bins show qualitatively the same trends. 

We find no strong correlation between environment and 
spin - it is possible to find low spin haloes in very dense en- 
vironments, and high spin haloes in sparse environments. 
However, we do find a small excess probability that haloes 
with lower spins are in less dense environments, as evi- 
de nced by the histog r ams in Figure (5] This is in contrast 
to lAvila-Reese et al.r (|2005l ). who report that galaxy scale 
haloes (masses < 5 x 10^^ M0/h) in clusters have lower spin 
than isolated haloes in the field. Interestingly, they find that 
the trend they report also disappears as they go back to 
2 = 1. This may imply a trend with redshift where at high 
redshift, high spin haloes are in denser environments, and at 
very low redshift, high spin haloes are in less dense environ- 
ments. However, a strict comparison of these two numerical 
studies cannot be made as we have significantly more haloes, 
probe significantly higher redshift sli c es an d significantly 
lower halo masses. lAvila-Reese et al] l|2005h only use the 
over-density criterion to quantify the environment, whereas 
our findings hold across three different measures of environ- 
ment. We note that the trends reported here have shallow 
slopes, similar to the lMaccio et al.1 (|2007l ) work. We conclude 
that the dependence of spin on environment is weak for the 
halo masses and epochs studied here. 



6 STRUCTURAL PROPERTIES 
6.1 Concentration 

In determining the effect of halo spin on halo structural 
properties, we first turn to the mass distribution within the 
halo. We fit the NFW profile, 



p(r)i. 



Po 



(4) 



r/va (1 + r/rs)^ 

(|Navarro. Frenk. fc Whitelll996l ■ ^1997^ ■ to each halo through 
a least squares fit to the radially averaged density profile 
calculated using all particles assigned to the halo by our 
group finding algorithm. In doing the least squares fit, we 
seek to minimize x^, which is defined as 



^ iVbi„= 

= - [log p - log PNFW] 

u ^ — ' 



(5) 



and where v is the number of degrees of freedom. We then 
use the concentration parameter, defined as Cns ~ Rns/rs, 
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Table 2. Statistical comparison of our LoRes and MedRes runs. We present the mean and standard deviation of seven quantities: the 
offset between the centre of mass (in physical kpc); the fractional difference (defined as \Xm — Xi\IXi for property X) in mass, spin, 
kinetic energy, potential energy, and total angular momentum; and the cosine of the angle between the angular momentum vectors in the 
two runs. We find that the scatter decreases with increasing resolution for all halo properties, and that the angular momentum vector is 
more robust at large particle numbers. 



Property 


N < 1000 


1000 <N < 3160 


3160 < Af < 10000 


N > 10000 




Mean 


(J 


Mean 


fj 


Mean 


fj 


Mean 


a 


Ar 
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0.162 


0.185 
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0.158 


0.173 


0.171 


0.198 


AM/M; 


0.141 


0.176 


0.107 


0.114 


0.072 


0.119 


0.046 


0.066 


AA/A, 


0.563 


0.873 


0.422 


0.692 


0.257 


0.409 


0.254 


0.532 


Ai?K/(i?K)i 


0.309 


0.483 


0.168 


0.166 


0.102 
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0.097 


0.226 




0.286 
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0.094 


0.158 


0.102 


0.230 
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1.24 


0.377 


1.85 


0.248 


0.584 


cos{e) 
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0.475 


0.745 


0.388 


0.857 


0.300 


0.916 


0.179 
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Figure 5. Two measures of environment versus spin parameter: the distance to the 7"* nearest neighbor {Dj, left column), and the 
over-density within lO-Rirg (right column). The solid curves show the mean spin, and the error bars reflect Poisson noise. The top two 
rows show two different mass bins. The bottom row plots histograms of the 10^ halo sample, binned by spin such that one third of 
the haloes lie in each of t he hig h and low spin bins. Wc find only weak trends in each mass bin, but they do confirm our main finding 
from lDavis fc Nataraiaiil l|2009l ') that higher spin haloes are more likely found in denser environments as they have smaller values of D7 
and reside in over-dense regions. 



to characterize the density profile. In Figure |6] (left column) 
we show the measured concentrations versus spin parame- 
ters binned by mass. There is only a small trend with spin, 
implying that the density structure has only a weak depen- 
dence on the angular momentum properties of the halo. In 
the bottom left panel of Figure [S] we show the histogram of 
concentration values for haloes in the 10*' Mq mass bin. We 



report little to no dependence of C178 on spin in these mass 
ranges. We also look at the relationship betw een halo mass 
and Ci7 8 . We find that the fitting function in iBuUock et al.l 
extrapolated down to these masses and up to these 
high redshifts, given by Cvir « 9m-° "/(1 + z), where 
fi — M/Af* and is the typical collapsing mass at redshift 
z, is not a good fit, and predicts values for the concentra- 
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tion that are too small compared to our estimates from the 
simulation. This disagreement is not surprising as Bullock 
et al. studied higher mass haloes at lower redshifts {z < 5) 
that have had a significant amount of time to merge and 
virialize. 



6.2 Circular Velocity 

Observed rotation curves are used to obtain mass estimates 
of gala^xies and clusters. The measured velocity of stars and 
gas reflects the combined gravitational potential of both the 
baryons and dark matter in the galaxy. To infer the to- 
tal mass from an observationally determined rotation curve 
requires understanding how both types of matter are dis - 
tributed spatially in galaxies. IPersic. Salucci. fc Steil (|l996l ') 
report that for Sb-Im spirals, rotation curves can be repre- 
sented by a universal function which is the sum in quadra- 
ture of two velocity curves: one from the disc and one from 
the dark matter ha lo. This model is developed further in 
ISalucci et all (|2007l '). where the dark matter velocity com- 
ponent is strictly a function of the virial mass of the halo. 
If, however, the velocity curve of the dark matter halo is de- 
pendent on a second parameter, such as A, we would expect 
to find a systematic error in a mass estimate for the galaxy 
derived from the observed rotation curve. 

A relationship between A and a characteristic circular 
velocity can be ex pected if we consider an alternative spin 
parameter given in lBullock et all (|2001B| '): A' = J/V2MVR, 
where J is the angular momentum inside some radius, R, M 
is the mass inside R and V is the circular velocity at R. For 
the case of a truncated, singular isothermal sphere. A' = A 
at the virial radius. Thus it is expected that a relationship 
between the A used in our work and a characteristic circular 
velocity should exist. 

We find that the peak velocity, Knax, in the meas ured 
circular velocity curve, defined as Vc — ^JGM{< r)/r, sys- 
tematically depends on spin. To measure this curve, we sort 
all particles in the halo by their distance to the halo centre. 
From this, we can read off the circular velocity curve at the 
radius of each particle, and we define Vmax as the maximum 
value of this curve. Figure [6] (right column) shows that for 
all mass bins, there is a systematic offset in the peak veloc- 
ity as a function of halo spin: high spin haloes have higher 
Vmax than their low spin counterparts. The bottom panel 
shows two histograms of Vmax, binned by spin. The peak of 
the high and low spin curves is offset by 20% in the 10^ M0 
mass bin, and 10% in the 10^ M0 haloes (not plotted). 

From a given Vmax one can infer a halo mass using 
the definition of circular velocity such that M(< r) — 
Vnax^'/G. Alternatively, one can use published relationships 
between Vmax and Mvir such as the power l a w Vn ax oc 
M^/ " where a = 0.31 ± 0.08 fro m 'Shaw et al.' (2006') (see 
also Kravtsov, Gnedin. fc Klvpin 2004; Hayashi ct al. 200^1; 



iBuUock et al.ri2001Al ). Either way, when two halos of the 
same (unknown) mass have different spins, there is a bias 
in the resulting Vmax, which will lead to a bias in the mass 
estimated from the circular velocity. Therefore two ob- 
served galaxies with differing values of Vmax could 
in fact inhabit dark matter haloes of the same mass 
but with different spin parameters. Alternately, two 
galaxies with the same value of Vmax could have dif- 
ferent masses due to their different spin values. Wc 



can estimate the extra uncertainty in the halo mass due to 
the bias induced by A. Using the circular velocity definition, 
we can translate a difference in velocity to a difference in 
mass: 2A{v)/v = A(M)/M. Therefore the estimated total 
mass from an observationally measured value of Vnax can be 
off by 20-40% arising due to the unknown value of the spin 
parameter of the dark matter halo that hosts the galaxy. 

6.3 Halo Shape and Spin 

Another well explored correlation at low redshift is between 
the shapes of dark matter haloes and the spin parameter. 
Here, we investigate the corresponding correlations for these 
high redshift haloes. The existence of a correlation between 
halo shape and spin is likely to have important consequences 
for the formation of pre-galactic disks in these haloes. 

We calculate the sphericity, s = c/a, and triaxiality, 
T = (a^ — b'^)/{a'^ — c^), from the eigenvalues (a > 6 > c) of 
the normalized moment of inertia tensor, 

I.J 



EXiXj 



where x is the distance to the halo center. The ten- 
sor is calculated using only the particles assigned to the 
halo, not from all particles within the i?i78 as is often 
done. This allows for consistency with our measurements 
of A. We used the normalized tensor so that we do not 
weight particles on the outskirts of a halo stronger. This 
helps to prevent the halo shape being dominated by resid- 
ual tidal features from a recent merger or other irreg- 
ularities found only in the outer regions which are not 
present in th e halo interior. We caution th at both the nor - 
malized fe.g. lAvila-Reese et al.l[2 005: AUg ood et aP I2OO6I) 
and unno rmalizcd (e.g. Jans-CondcU & Hcrnauist 200lJ; 



Shaw et al. 2006: Maccio^ Duttori, & van den Bosch 2008 



Faltenbacher fc Whit3l2O10l l tensor are used in the litera- 



ture, and so caution must be used when comparing shape 
distributions. We show in Figure [7| histograms of s and T 
when calculated using both a normalized and unnormalized 
moment of inertia tensor. We find that many more halos 
are significantly aspherical, and halos tend to be less pro- 
late when the shape is found usi ng the unnormalized tensor. 
This is contrary to the study in (|Allgood et al.|[2006l ) , which 
reported no systematic offset in s between the two meth- 
ods. We note, however, that they only used particles within 
O.Siivir, rather than all particles assigned to the halo as we 
do. Therefore, our unnormalized method will have a stronger 
bias because we include distant particles in the calculation 
of the halo shape. In what follows, we use the normalized 
tensor to calculate the shape of the haloes. 

We show in Figure |8] the relations between s (left col- 
umn), T (right column), and A, as well as the histograms of 
s and T binned by spin and mass. For higher spin haloes, 
the trend is for haloes to become less spherical, and more 
prolate. This is in agreement with studies at low redshift 
and higher rnasses (IMaccio. Dutton. fc van den Boschll2008l: 
'Bett et all l2007l : lAUgood et all I2OO6I : lAvila-Reese et all 
2005). 

We also plot the two point correlation function, ^(r) of 
the haloes in Figure [9l binning the haloes by their values 
of s and T. The correlation function represents the excess 
probability of finding a halo at a distance r when compared 
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Figure 6. Concentration (left column) and Vinax/ < V > (right column) versus spin for haloes binned by mass. Vmax is the maximum 
value of the circular velocity curve, given by Vc = yjGM{< r)/r, and < y > is the mean value of Knax for all haloes in the given mass 
bin. The bottom row shows histograms of the 10® Mq haloes, binned again by spin. While the concentration does not vary much with 
spin, the peak of the rotational curve does vary with spin. This implies that deriving masses from observed velocity curves will have an 
additional systematic error due to the spin of the host dark matter halo. There is a difference of 20% in the peak value for 10® Mq haloes 
and a 10% spread for the 10^ Mq haloes with spin parameter. 



5000 ^ 




Figure 7. Histograms of sphericity, s (left), and triaxiality, T (right), for all haloes in the WM5 run. We show the values of s and T 
found using a normalized (solid black curve) and an unnormalized (dashed blue curve) moment of inertia tensor. We find that using the 
unnormalized tensor gives significantly more aspherical haloes and fewer prolate haloes than using the normalized tensor. 
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Figure 8. Sphericity, s (left panels), and triaxiality, T (right panels), versus spin parameter for haloes binned by mass. High spin haloes 
are less circular than their low sin counterparts. High spin haloes are also more likely to be prolate. 



to a random distribut ion of haloes. We use the same method 
as our previous work (jPavis fc Nataraiarj [2009') to calculate 
^(r). We calculate ^(r) for haloes in four different bins, cor- 
responding to large and small values of s (left column of 
Figure O and T (right column) . The cuts were chosen so 
that one third of the haloes lie in each bin. We also separate 
out haloes by mass so that we eliminate any mass effect on 
the value of ^(r). In the top row of Figure|5]we show ^(r) for 
haloes with M = lO'^ Mq and in the bottom row we show 
10'' M0 haloes. 

We find that nearly spherical haloes (values of s close to 
1) are more clustered than the more aspherical haloes. This 
trend is stronger for haloes in the 10^ M0 mass range than 
in the 10^ Mq mass range. In the higher mass bin, we find an 
increase of 50% in the correlation of nearly spherical haloes 
when compared to the extremely aspherical halo sample. In 
the lower mass range, we report a smaller increase of only 
20%. Thus we conclude that our r esults at high redshift 
follow trends f ound at lower redshift fFaltenba cher fc White! 
l2010l : lBett eT al. 2007; Avila-Rccsc ct al. 2005). 

Unlike the sphericity parameter, when we separate 
haloes into high and low values of T, we see little change in 
^(r). T his is somewhat surprising, as Fal tenbacher fc White! 

find at low redshift an offset in ^(r) when binning 
haloes by T. Finally, we note that the clustering strength is 
weaker if an unnormalized moment of inertia tensor is used 
to calculate s and T of the haloes. 



In Figu re IIOI we show the bias parameter, b = 
\J ^MM / ^HH , as a function of the halo peak height 
v(M, z) = &cl{( y{M)D{z)), where D{z) is the growth func- 
tion (|Mo fc White 2002 ). This allows us to compare results 
from two redshifts (2 = 10 and 2 = 6) and from 5 mass 
bins at each redshift. We split our haloes according to four 
properties: A, s, T, and Cits- We report that the spin pa- 
rameter has the largest effect on the bias, while the triax- 
iality has the least. This implies that angular momentum 
has the strongest dependence on environment of all these 
considered variables and will be the source of the largest 
systematics due to differing environments between haloes. 
Therefore, baryonic properties that depend on the angular 
momentum of the host halo should have systematic offsets 
due to their local environment. This may play a large role in 
the formation of the earliest galactic disks, as well as affect 
semi-analytic models which relate the ba r yonic spin to the 
dark matter halo spin ICroton et al.ll2006l : iBenson fc Bowen 
!2010l , e.g.), and estimates of t he dark matter halo spin 
from the observed baryonic disk (Cervantes-Sodi et al.!!2008! : 



ICervantes-Sodi. Hernandez fc Park 2010l ). 



7 DISCUSSION AND CONCLUSIONS 

Our key findings can be summarized as follows: 

• We have measured the spin, concentration, circular ve- 
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Figure 9. Correlation function for haloes binned by their sphericity (left), triaxiality (right) and mass. We find that nearly spherical 
haloes are more clustered than aspherical haloes, by 50% in the 10^ Mq bin and by 20% in the 10*' Mq bin. We find no difference in 
clustering strength when separating haloes based on T. 



locity, sphericity, and triaxiality parameter for a statistically 
large sample of dark matter haloes at high redshift {z > 6). 

• High spin haloes at high redshift are 25% more clustered 
than their low spin counterparts at a given mass, and are 
more likely to be found in high density environments. 

• High spin haloes (with masses < 10^ M©) have smaller 
maximum circular velocities than low spin haloes, leading 
to errors up to 40% in the derived enclosed mass. 

• High spin haloes at high redshift are more likely to be 
aspherical and prolate, similar to findings at low redshift. 

• Nearly spherical haloes are up to 50% more clustered 
than extremely aspherical haloes, while there appears to be 
no difference in the clustering strength based on the triaxi- 
ality of the haloes. 

Our findings have an impact in two general areas: the 
role of angular momentum in halo structure and formation, 
and the role of assembly bias at high redshift. Our find- 
ings show that angular momentum has a measurable cor- 
relation with structural properties, including the concentra- 
tion, sphericity, and triaxiality. Also, haloes with higher spin 
are preferentially found in higher density environments. The 
finding that halo spin correlates with local environment at 
high redshift is important to the understanding of the evolv- 
ing properties of the baryonic component of dark matter 
haloes. A correlation between spin and baryonic properties, 
such as formation time, disk rotational speed, or disk size. 



would be the specific consequences of the correlation with 
the environment. These correlations are likely to be signif- 
icantly stronger at high redshift, before too many mergers 
have happened which could destroy any correspondence be- 
tween the dark matter spin and the baryonic structure. 

The role of assembly bias in halo evolution has been 
discussed at low redshifts. Our work extends the study of 
assembly bias to high redshifts, when the first galaxies form. 
We see similar results to those at low redshift when studying 
the dependence of clustering on spin and sphericity of dark 
matter haloes. However, our results differ when looking at 
the triaxiality of haloes. We find that there is no difference 
in clustering strength between prolate and oblate haloes. 

Our findings are of particular importan ce now that 
galaxies are being found at these high redshifts (|Oesch et al.l 
2010). The clustering of galaxies has been use d to infer the 
masses of their host dark matter haloes (e.g., lOuadri et al.l 
120081 ). However, we find that properties other than mass 
- in particular halo spin - affect the measured correlation 
function. This additional parameter will induce errors in 
mass estimates for dark matter haloes inferred purely from 
clustering measure ments. This fits in with the results of 
lOuadri et al.l (|2008l ). who suggest that mass is likely not the 
only parameter that drives the interaction of haloes with 
their large-scale environment. 
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Figure 10. Bias as a function of peak overdensity, u{M,z). We calculate the bias at two redshifts: z = 6 (filled symbols) and z = 10 
(open symbols). The top left panel bins haloes by their spin parameter, the top right by their sphericity, s, the bottom left by the 
triaxiality, T, and the bottom right by the concentration, Cng. We see a strong offset in bias due to spin and concentration, a weaker 
one due to s, and no offset due to T. 



In addition to mass measurements, assembly bias will 
play an important role in feedback at these high redshifts. 
At the highest redshifts, simulations show that Population 
III stars have a large impact on t heir environment due to ra- 
diative and supernova feedback (jjohnson. Grief, fc Bromml 
l2007l : iGrief et ahlboOTl : IWhalen et al.ll2008 a"b'). One conse- 
quence of our findings is that if angular momentum affects 
the formation and evolution of these Pop III stars, their 
feedback effects will show an environmental bias. Thus, the 
distribution of metals and reionization will be more clus- 
tered than otherwise expected. We also expect that because 
these stars are the first baryonic objects to collapse, their 
properties can be expected to have a stronger relationship 
to their host dark matter halo than galaxies today, which 
have undergone multiple mergers. We intend to pursue the 
consequences of our findings on baryonic results in future 
work. Our results suggest that the angular momentum prop- 
erties of dark matter haloes likely have consequences for the 
properties of the first stars and galaxies hosted by them. 
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